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Abstract. We show that generalised extreme value statistics —the statistics of 
the k th largest value among a large set of random variables- can be mapped onto 
a problem of random sums. This allows us to identify classes of non-identical and 
(generally) correlated random variables with a sum distributed according to one 
of the three (fc-dependent) asymptotic distributions of extreme value statistics, 
namely the Gumbel, Frechet and Weibull distributions. These classes, as well as 
the limit distributions, are naturally extended to real values of k, thus providing 
a clear interpretation to the onset of Gumbel distributions with non-integer index 
k in the statistics of global observables. This is one of the very few known 
generalisations of the central limit theorem to non-independent random variables. 
Finally, in the context of a simple physical model, we relate the index k to the 
ratio of the correlation length to the system size, which remains finite in strongly 
correlated systems. 
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1. Introduction 

One of the cornerstones of probability theory is the celebrated central limit theorem, 
stating that under general assumptions, the distribution of the sum of independent 
random variables converges, once suitably rescaled, to a Gaussian distribution. This 
theorem provides one of the foundations of statistical thermodynamics and, from a 
more practical point of view, legitimates the use of Gaussian distributions to describe 
fluctuations appearing in experimental or numerical data. 

However, this general theorem breaks down in several situations of physical interest, 
when the assumptions made to derive it arc no longer fulfilled. For instance, if the 
random variables contributing to the sum have an infinite variance, the distribution 
of the sum is no longer Gaussian, but becomes a Levy distribution pp. This has 
deep physical consequences, for instance in the context of laser cooling [5] or in that 
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of glasses j3] 0] . Though the central limit theorem applies to a class of correlated 
variables, the martingale differences [3] E], it is generally inapplicable to sums of 
strongly correlated, or alternatively, strongly non-identical random variables. In 
this case, no general mathematical theory is available, but experimental as well as 
numerical distributions of the fluctuations of global observables in correlated systems 
often present an asymmetric shape, with an exponential tail on one side and a rapid 
fall-off on the other side. Such distributions arc usually well described by the so-called 
generalised Gumbcl distribution with a real index a. The latter is one of the limit 
distributions appearing, for integer values a = fc, as the distribution of the k th extremal 
value of a set of independent and identically distributed random variables Q |S| ■ For 
instance, Bramwell and co-workers argued that the generalised Gumbel distribution 
with a ss 1.5 reasonably describes the large scale fluctuations in many correlated 
systems jH]- Since then, the generalised Gumbel distribution has been reported in 
various contexts ^3^2 El E3 EH Yet, the interpretation of such a distribution 
in the context of extreme value statistics (EVS) is far from obvious. This issue led 
to a debate around the possible existence of a (hidden) extreme process which might 
dominate the statistics of sums of correlated variables in physically relevant situations 
|15l llfil rm I18j . Still, no evidence for such a process has been found yet. And even if 
one accepts that some extremal process is at play, a conceptual difficulty remains, as 
speaking of the fc th largest value for non-integer k simply does not make sense. 
In this article we propose another interpretation of the usual EVS in terms of statistics 
of sums of correlated variables. Within this framework the extension to real positive 
values of k is straightforward and avoids any logical problem. Usual EVS then appears 
as a particular case of a more general problem of limit theorem for sums of correlated 
variables belonging to a given class. 

The article is organised as follows. In Section 2 we present the link between usual EVS 
and sums of correlated variables, and propose a more general problem of random sums, 
which is studied is the next sections. In Section 3, we present in detail the simple 
problem of sum -already presented shortly by one of us in a previous publication 
|19|- associated with the EVS of exponential random variables. Section 4 deals with 
the extension to the more general problem of random sums defined in Section 2, 
and it is shown that the asymptotic distribution of the sum, is either a (generalised) 
Gumbel, Frechet or Weibull distribution. Finally, we discuss in Section 5 some physical 
interpretations of our results. 

2. Equivalence between extreme value statistics and sums of correlated 
variables 

2.1. Standard extreme value statistics 

The asymptotic theory of EVS, which found applications for instance in physics |2(Jj . 
hydrology J2J| ; seismology [22] an d finance , is also part of the important results 
of probability theory. Consider a random variable X, distributed according to the 
probability density P, and let jEjv be a set of -/V realisations of this random variable. 
Now instead of considering the sum of these N realisations as in the central limit 
theorem, one is interested in the study of the distribution of the fc th maximal value 
Zk of Xjv- The major result of EVS is that, in the limit N — > oo, the probability 
density of does not depend on the details of the original distribution P: There are 
only three different types of distributions, depending on the asymptotic behaviour of 
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P [3IH1- In particular, in the case where P decreases faster than any power laws at 
large x, the variable z k is distributed according to the Gumbel distribution 



^ = *'(*), VH = 



-k9 k (z k + v k ) - fc e -M*»+>*) 
J-[lnA -*(*)], 



(1) 



where ^(x) = ^ lnT(a;) is the digamma function. Note that from the point of view 
of EVS, k G N* by construction. Yet, if one forgets about EVS, the expression is 
also formally valid for k = a real and positive. 

2.2. Reformulation as a problem of sum 

Let us consider again the set Xn defined here above. We then introduce the ordered 
set of random variables, z k = x (T (fc), where a is the ordering permutation such that 
Z\ > %2--- > Zn- One can now define the increments of z with the following relations: 
u k = z k - z k+u Vl<fc<iV-l, 

UN = Z N . 

Then by definition, the following identity holds 



Max 



(") 



N 



Uk, 



(2) 



where Max'"' is the n th largest value of the set Xn , showing that an extreme value 
problem may be reformulated as a problem of sum of random variables U k . 
Besides, as these new random variables have been obtained through an ordering 
process, they are a priori non-independent and non-idcntically distributed. This is 
confirmed by the computation of the joint probability J k N{u k , ...,Un) which is given 

by 



Jk,N{uk,:;U N ) = N\ dz N P(z N ) / dz N -iP(z N -i) 



dziP(zi) 



JV-l 



x S(u N ™ z N ) \\ 8(u n ~ z n + z n+ i), 

n—k 

\P (J2 Ui ) / dzk-iP(zk-i)... d Zl P( Zl ). 

n=k \i=n J \ Ju k +...+u N Jz2 

F(u) k -\ 



= TV! 



By recurrence it is easy to demonstrate the following relation 

1 



p oo poo 

/ dzk-iP{zk-i)- / dziP(zi) 

J U J Z2 

with F(u) = dz P(z), leading to 

/ AT 



(fc-1)! 



fc-i 



JV 



N 



Jk,N(u k , ...,U N ) 



(fc-1) 



The probability J k: N is actually a probability at N 1 = N + 1 — k points, and a shift 
of indices allows one to write it in the final form: 

fc-i 

(TV' + fc-l)! 



N' 



N 



JV 



Jn'{u\ 



.,Un> 



(fc-1)! 



(3) 
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In the general case, this expression docs not factorise, i.e. one cannot find N' functions 
7r„ such that 

N' 

Jn'(u%, ...,Utf') = Y\_ KniUn), 
n = l 

and the random variables U n are thus non-independent. 
2.3. Extension of the joint probability 

Up to now, we only reformulated the usual extreme value statistics as a problem of 
sum, by some trivial manipulations. Basically, wc have two equivalent problems with 
the same asymptotic distribution, e.g., the standard Gumbcl distribution if P is in the 
Gumbel class. Let us now point out that Eq. is a particular case of the following 
joint distribution 

j N ( Ul ,...,u N ) = ^n ^(E u ") Il p (j2 u ) ( 4 ) 

\n=l / J n—1 \i—n / 

where O(-F) is an (arbitrary) positive function of F, and where Zn is given by 



Z N = / dv£l(v){l-v)"-\ (5) 



o 



Eq. © is recovered by choosing Q(F) = F . Note that by extending the definition 
of the joint probability, we lost the equivalence between the two original problems. In 
other words, we are generalising the problem of statistics of sums and not the extreme 
value one. In the rest of this article, we study the limit distribution of a sum of 
correlated random variables satisfying the generalised joint probability Q. 
Let us note finally that one can easily generate from the joint probability (0J other joint 
probabilities that lead to the same asymptotic distribution for the sum, by summing 
over a set S of permutations over [1, N] 



1 



J%(ui, ...,Un) = T77c\ J n(Uct(1), Ua(N)), 

^ ' aes 

J\f(S) being the cardinal of S. This may allow in particular some symmetry properties 
between the variables to be restored. Indeed, starting for instance from independent 
and non-idcntically distributed random variables, this procedure leads to a set of non- 
independent and identically distributed random variables with the same statistics for 
the sum. 



3. The exponential case 

Though in general the joint probability docs not factorise, there is a particular case 
where it does. For pedagogical purposes, wc discuss this case separately, before dealing 
with the general case in Section^] To allow the joint probability J/v defined by I0J to 
factorise, the function P has to satisfy the following property: 

V(a?,y), P(x + y) = P(x)P(y). 

That is to say that P must be exponential: P(x) = ne~ KX , with k > 0. As a result, one 
has F(x) = e~ KX , and the factorisation criterion on Jn translates into the following 
condition on Q 

V(F U F 2 ), nCFiJb) = n(F 1 )Q(F 2 ), 
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so that f2 must be a power law: fl(F) = F -^ 1 (the prcfactor may be set to 1 without 
loss of generality), with a > to ensure the convergence of the integral defining Zn 
given in Eq. Ipjl. In such a case, one finds 

Zn T(N + a)' (6) 
and the joint probability could be written in the factorised form Q, with 

7r n K) = (n + a-l) Ke -(" +a - 1 ) K "". (7) 

Therefore, for a = k integer, if one sums independent and non-idcntically distributed 
random variables U n , n = 1, N, obeying 0, the distribution of the sum converges 
to a Gumbel distribution Gk in the limit N — > oo. Note that the particular case k = 1 
has been previously studied by Antal et al in the context of 1// noise [2I|. In the 
following of this section, we shall establish that for a real, the distribution of the sum 
of the U n 's converges towards the generalised Gumbel distribution G a for N — > oo, as 
announced in a previous publication |19| . 

To this aim, let us define the random sum Sn = J2n=i Un, where U n is distributed 
according to 0, with a real and positive. The distribution of Sn is denoted by Tjv- 
As the Un's are independent, the Fourier transform of Tn is simply given by 

N N 



n(n + a — 1) / 



>N 



n—1 n— 

The first two moments of the distribution read 

N N 

(s N ) = x>„> = -x;— — r, (9) 

^— 4 k Z — ' n + a — 1 

n— 1 n—1 

iV AT 

53 Var(tfn) = -jE " ( 10 ) 

n=l n=i v ' 

In order to get a well-defined distribution in the limit N — > oo, let us introduce the 
reduced variable \i by 

s - (S N ) .., ,. 
/i = , with cr = lim (Tat. 

a N-^oo 

Note that the fact that a < oo breaks the Lindeberg's condition, allowing an eventual 
breakdown of the central limit theorem |2S| . The distribution of /x, $at, is then given 
by &n(h) = o~Tn(o~h + (Sn))- The Fourier transform of $ = limTv^oo 'I'iv reads, 



using Ij8l9jl 

[$„,] (w) = lim ^ [$at] (w) = TT f 1 + 7 T tt) exp f -) 

n^oo „- 1 V cr«(n + a — 1) / \aK(n + a — l)y 

The last part of the previous expression has to be compared with the Fourier transform 
given in Appendix A, leading to 

$oo (M) = <3 (/i). 

This result shows that it is possible to obtain quite directly the generalised Gumbel 
distribution with a real index a, not from an extremal process, but from a sum of 
independent non-identically distributed random variables. As already pointed out, the 
fact that the random variables U n are uncorrelated is a specificity of the exponential 
distribution. A natural question is then to know whether this result survives for a 
more general class of distributions P and functions 0. We address this question in 
the next section. 
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In the general case, the random variables U„ are non-independent, so wc have to 
deal with the joint probability (@J. Accordingly, the route to the asymptotic limit 
distribution will be quite different from the exponential case. 



4-.1. Distribution of the sum for finite N 

Consider a set of realisations {u n } of N (correlated) random variables U n , with the 
joint probability |0J. We then define as above the random variable Sn = Y^ n =i 
and let T^r be the probability density of Sn- Then Tjy is given by 



T N (s) = 
Inserting Q}, one obtains 



,-oo / n n 

/ diijv.-dui Jn(ui, un) 6 I s — u n 

J ° V n=l , 



= ^l P (s)a(F(s))lN(s), 



with 
In(s) 



du N P(u N ) du N ^iP(u N + u N -i). 



r°° ( N \ 

J o \ „=i / 



To evaluate In, let us start by integrating over u±, using 



J du^ (s - ^2 U nj = 9 ^ S _ X] ' 



where is the Heaviside distribution. This changes the upper bound of the integral 



over U2 by u\ 



max 



(o,*-Et 



_2 "n 



Then the integration over U2 leads to 



N 



\n— 3 



du 2 P(u2) 

By recurrence it is then possible to show that 
In(s) 



F E Un - 



iV 



e 



1 / N^- 1 

r(7V) V v o 

finally yielding the following expression for T n ■ 

1 / \N-1 

T n(s) = -^P(s)a(F(s)) (1- F(s)) . 

In the following sections, we assume that £l(F) behaves asymptotically as a power 
law £l(F) ~ iloF"" 1 when F — > (a > 0). Under this assumption, wc deduce 
from Eq. the different limit distributions associated with the different classes of 
asymptotic behaviours of P at large x. 



(11) 
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4-2. The Gumbel class 

In this section we focus on the case where P is in the Gumbel class, that is, P(x) 
decays faster than any power law at large x. Note that the exponential case studied 
above precisely belongs to this class. Our aim is to show that, after a suitable rescaling 
of the variable, the limit distribution obtained from Ijlljl is the generalised Gumbel 
distribution G a , where a characterises the asymptotic behaviour of Q(F) for F — > 0. 

To that purpose, we define by F(s^) = a/N . If a is an integer, this is nothing 
but the typical value of the a th largest value of s in a sample of size N. As P is 
unbounded we have 

lim s* N = +oo. 

N— too 

Let us introduce g(s) = — hiF(s) and, assuming g'(sisr) ^ 0, define the rescaled 
variable v by 

s = s * N + 7k~Y (12) 

For large N, one can perform a series expansion of g around s 

^ 1 „ 



* . 
jv- 



g(s) = g(s* N ) +v + 



~ x n\ g'(s* N Y 



For P in the Gumbel class, g^ n ' (s^) / g' (s* N ) n is bounded as a function of n so that 
the series converges. In addition, one has 

lim 4^=0, Vn>2, 

so that g(s) may be written as 

g(s) ~ g(s* N ) + v + Sn(v), with lim En(v) = 0. (13) 

N^oc 

Given that P(s) = g'(s) F(s), one gets using Eqs. ifTTJl and (|TS|> 

# » w " i4) TK(,) 

= ihMj FmiF{s)) ( 1 - F(s) Y~ 1 - 

where s is given by Eq. (|12(l . For P in the Gumbel class, it can be checked that, for 
fixed v 

N->oo g'( S * N ) 

Besides, F(s* N + v/g'(s* N )) -> when N — > oo, so that one can use the small F 
expansion of Cl(F). Altogether, one finds 

a 



fa\ a _„,,_ nFKjhA \ _ a_ v _ En{v) -\ r 
N 



—av—aett(v) 



1 ' z N \nJ 

Using a simple change of variable in Eq. , one can show that 
J lm — o — = r W 

N — >oo SZq 

It is then straightforward to take the asymptotic limit N — > oo, leading to 

a a 

*oo(«) = fTTT ex P [~ av ~ a e~ v ] ■ 

r(o) 
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Gumbel G k , ke N* 



Gumbel G a . a t ] 



Random sum with J* 



■ Random sum with J„ 



Figure 1. Generalisation of EVS via sums of non-independent variables. 



In order to recover the usual expression for the generalised Gumbcl distribution, one 
simply needs to introduce the reduced variable 

v — (v) 
fi = , 

with, ^ being the digamma function, 

(v) = Lna-tf(a), o% = tf'(a). 

The variable fi is then distributed according to a generalised Gumbcl distribution Q . 
To sum up, if one considers the sum Sn of iV ^> 1 random variables linked by the joint 
probability then the asymptotic distribution of the reduced variable ^ defined by 

s N - {S N ) 



(TJV 



with 



(S N ) = 



hu 



*(o) 



'JV 



0~N 



g'( s * N ) ' 9'{s* N ) ' 

is the generalised Gumbel distribution QJ. More generally, the approach developed in 
this paper to relate EVS with sums of non-independent random variables, and then 
generalise the problem of sum, may be summarised as shown on Fig.^on the example 
of the Gumbel class. 



4-3. Extension to the Frechet and Weibull classes 

Let us now consider the cases where P belongs either to the Frechet class, that is P 
has a power law tail 

x — >oo X ' ' 

with P > a constant, or to the Weibull class, that is P has an upper bound A and 
behaves as a power law in the vicinity of A 

\fx > A, P(x) = 

P(x) ~ P (A-xf-\ /3>0. (14) 
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with P an arbitrary positive prefactor. A calculation similar to the one used in the 
Gumbel case allows the asymptotic distributions to be determined. Considering first 
the Frechet class, we define through F(s* N ) = a/N as above, and introduce the 
scaling variable v = s/s* N . This gives for the distribution $> l N (v) 



v 



N a Z N v 1+a7 i 



1 - 



N-l 



Taking the limit N — > oo, the asymptotic distribution reads 



V 



-a/v" 



T(a) v 1+a7 i 

One can rescale the variable v with respect either to the mean value or to the standard 
deviation a. In this latter case, one finds 



,l/v 



Introducing (i 



T(a) 
■ vja we get 



r(a)r (a--)-r(a-- 



1/2 



1 



T(a) fi 
r(a)" 



exp(-A^ r ') 



r(a)r [a- - ) -V la- - 



i) 



V 



1 -»7/2 



Calculations for the Wcibull class are very similar, given that the scaling variable is 



(A- 



s)/{A 



Skipping the details, the asymptotic distribution reads 



f3a a 



a/3- 



After rcscaling to normalise the second moment 
a" 1 '? 



T(a) 



r(o)r 



I 1 



-,1/2 



one obtains 



P 

r(a) 



A 11 ^- 1 exp (-A/) , 



r(o)r 



a+ - - 



Thus here again, the distributions appearing in EVS are generalised into distributions 
of sums, including a real parameter a coming from the asymptotic behaviour of the 
function fl. 

Note that the generalised Frechet and Weibull distributions found above may also 
arise from a different statistical problem. As already noticed by Gumbel \7\ in the 
context of EVS, the Frechet and Weibull distributions may be related to the Gumbel 
distribution through a simple change of variable. If a random variable X is distributed 
according to a Gumbel distribution G a , then the distribution of the variable Y = e xx 
(A > 0) is a generalised Frechet distribution. As X is defined as a sum of correlated 
random variables associated with a joint probability J3J with P in the Gumbel class, 
Y appears as the product of variables with P in the Frechet class. Similarly, the 
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generalised Wcibull distribution may also be obtained from the Gumbel distribution 
by the change of variable Y = e~ xx . 

In summary, one sees that the reformulation of standard EVS by means of the joint 
probability allows a natural and straightforward definition of all the generalised 
extreme value statistics. As already stressed, this generalisation breaks the equivalence 
with the extreme value problem: there is no associated extreme process leading to the 
generalised extreme value distributions Gk, Fk or Wk with k real. At the level of the 
joint probability (@J, the equivalence with EVS only holds if Q(F) is a pure power- 
law of F, with an integer exponent k. Yet, concerning the asymptotic distribution, 
extreme value distributions are recovered as soon as Q(F) ~ fio F fe_1 for F — » (with 
k integer), that is, when is regular in F = 0. 

5. Physical interpretation 

Up to now we presented a mathematical result about sums of random variables, linked 
by the joint probability This joint probability leads to non-Gaussian distributions 
that may be interpreted as the result of correlations between those variables. All the 
informations about the correlation are included in Jn ■ From a physicist point of view, 
this is not completely satisfying: The degree of correlation is indeed usually quantified 
by the correlation length. At first sight it is disappointing that we can not extract 
such a quantity from |J3J. One has to realise however that up to now we only dealt 
with numbers, without giving any physical meaning to any of those numbers. To get 
some physical information from our result, we have to put some physics in it firstf, 
giving an interpretation of the u n , explaining also how those quantities are arranged 
in time or space, introducing then a notion of distance or time between the u n and the 
dimensionality of space. The results presented in the previous sections could therefore 
describe various physical situations. In this section we illustrate this idea by studying 
a simple model. 

Let us consider a one-dimensional lattice model with a continuous variable <f> x on each 
site x = 1, L. Although we do not specify the dynamics explicitly, we have in mind 
that the (j> x 's are strongly correlated. One can define the Fourier modes cf> q associated 
to <j) x through 



A natural global observable is the integrated power spectrum ("energy") E = ^2 X <fi x . 
From the Parseval theorem, E may also be expressed as a function of <j) q as E = 
J2 q l^gl 2 - We now assume that the squared amplitudes u n = |0 g | 2 , with q = 2irn/L, 
follow the statistics defined by Eq. (0J. For simplicity, we consider the simplest case 
where P(z) = Ke~ KZ and Q(F) = F 11 " 1 , although more complicated situations may 
be dealt with. This actually generalises the study of the 1// noise by Antal et al |24| . 
The power spectrum is given by 



f "Mathematician prepare abstract reasoning ready to be used, if you have a set of axioms about 
the real world. But the physicist has meaning to all his phrases." R. P. Feynman in 1261 . 



L 




(15) 
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The correlation function C(r), given by the inverse Fourier transform of \<fi q \ 2 , 

C{r)=T- 1 [(\^)]{r), (16) 

can be computed using 1)15(1 . leading to : 

C(r) oc — sin(mr) si(mr) — cos(mr) ci(mr), 

where si and ci are respectively the sine and cosine integral functions |27| . The 
correlation length is therefore defined as the typical length scale of C(r): 

€ = m" 1 = 



27r(o-l) 

Thus £ appears to be proportional to the system size, which was expected from the 
breaking of the central limit theorem \, The particular case of a 1// noise (a = 1) 
then corresponds to a highly correlated system, with £/L — > oo. For a sa 1.5, a value 
often reported in physical systems one gets £/L ~ 0.15: The correlation is weaker 
but £ still diverges with L [T8| . 

Altogether, this simple model may be thought of as a minimal model, that allows some 
generic properties of more complex correlated physical systems to be understood. 
Obviously not all correlated systems will exhibit generalised EVS: A well-known 
counter-example is the 2D Ising model at its critical temperature [2SI ■ Strictly 
speaking the generalised EVS can only be obtained if it exists P such as the joint 
probability is given by J3J. The number of observations of distributions close to a 
generalised Gumbo distribution suggests however that the expression Q is general 
enough to reasonably approximate the real joint probability in many situations. This 
would explain the ubiquity of the generalised Gumbel in correlated systems. 
More practically, the generalised extreme value distributions also appear as natural 
fitting functions for global fluctuations. If one is measuring fluctuations of some global 
quantities, it seems quite reasonable to fit them with an asymptotic distribution which 
could possibly take into account a violation of the hypothesis of the central limit 
theorem. The generalised EVS arc one of the few such distributions. Using the 
reduced variable fx, there is only one free parameter in the Gumbel class, a, which 
quantifies the deviation from the CLT. 



6. Conclusion 



In this article, we established that the so-called generalised extreme value distributions 
are the asymptotic distributions of random sums, for particular classes of random 
variables defined by Eq. Q), which do not satisfy the hypothesis underlying the central 
limit theorem. Interestingly, this is one of the very few known generalisations of 
the central limit theorem to non-independent random variables. In this framework, 
it becomes clear that it is vain to look for a hidden extreme process when one 
of the extreme value distributions is observed in a problem of global fluctuations. 
Therefore qualifying such distributions of generalised extreme value statistics is 
somehow misleading. The parameter a quantifies the dependence of the random 
variables, although further physical inputs arc needed to give a physical interpretation 
to this parameter. Within a simple model of independent Fourier modes, a is related 

| Dividing a system of linear size L into (essentially) independent subsystems with a size proportional 
to £, the number of subsystems remains finite when L — > oo if £ ~ L, so that the central limit theorem 
should not hold. 
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in a simple way to the ratio of the correlation length to the system size, a ratio that 
remains constant in the thermodynamic limit for strongly correlated systems. Besides, 
we believe that the classes of random variables defined by Eq. (0} may be regarded as 
reference classes in the context of random sums breaking the central limit theorem, 
due to their mathematically simple form. Along this line of thought, it is not so 
surprising that distributions close to the generalised Gumbel distribution arc so often 
observed in the large scale fluctuations of correlated systems, as for instance in the 
case of the XY model at low temperature. 

Appendix A. Fourier transform of the generalised Gumbel distribution 

In this Appendix we give an expression for the Fourier transform of generalised Gumbel 
distribution G a defined by 



r(a) 



-a9(x + v)- ae e{x+u) 



with 2 = ^'(a) and 6v = In a — ^(a). The Fourier transform of G a is defined by 
T [G a ] (w) = f + dx e- iwx G a (x). 



Letting u = acxp(— 6{x + v)), we obtain 

1 



r[G a ) H 



T(a) 

e ivu > 



a ' e 



du u 1 ' 



*/* r(a) 



Using the identity |2"7] 



ho— 1 — u 
e , 



r(i + «) = e -T*n r 



one finally obtains, with a = a — 1 



TT- 



exp 



[ (1 + a) iuJ / 

The identity HA.1|) leads to 

dlnr(l + a) _ 
da 



ioj s—^ f 1 1 

ft \ n n 4- 



-7 + 



n n + a 



L — ' \n n + a 



and therefore we obtain, using (|10|) 

a d 2 lnr(l + a) 
9 d^ 



oo ^ 

V _ 



(A.l) 



exp 



(n(0+a)) 



8(n+a) 



(A.2) 



a 2 . 



Furthermore, the relation (|A.2|I leads to 

°° / 1 1 
(„0 _ 7 ) - In (1 + <*) + £) 



n n + a 



[yd — 7) — In a + I 7 + 



dlnT 
da 
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So finally wc have 
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